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A HER X-l TURN-ON: USING THE PULSE PROFILE TO PROBE THE OUTER EDGE OF AN 

ACCRETION DISK 
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Abstract 

The X-ray binary pulsar Her X-l shows a wide vari- 
ety of long and short term variabilities in the X-ray light 
curve. The 35 d variability of the source is interpreted 
as the influence of a warped, inclined, and twisted ac- 
cretion disk periodically covering the line of sight to the 
central neutron star. In 1997 September we observed the 
"turn-on" of a 35 d cycle with the Rossi X-ray Timing Ex- 
plorer (RXTE) . Spectral analysis reveals that during early 
phases of the turn-on the overall spectrum is composed of 
X-rays scattered into the line of sight plus heavily ab- 
sorbed X-rays. This interpretation is consistent with the 
variation of the pulse profile observed at the same time. 
The overall shape of the pulse profile is not changing, but 
towards earlier phases of the turn-on the pulse signature is 
steadily "washed out" . This behavior can be understood as 
an influence of scattering and absorption due to the pres- 
ence of the accretion disk rim. Using a Monte Carlo code 
we simulate the influence of both processes on a time vari- 
able, beamed emission characteristic, similar to the pulse 
profile of Her X-l. By comparing the results of the simula- 
tion with the observed profiles we determine the amount 
of scattered radiation, absorbed radiation, and the size of 
the scattering region. 
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1. Introduction 

Her X-l/HZ Her is a X-ray binary system consisting of a 
neutron star Her X-l and its optical companion HZ Her. 
The system rotates with a 1.7 d orbital period around their 
common center of mass. The X-ray light curve of Her X-l 
shows strong variability on many time scales: e.g. 1.24 s 
pulsations originating in the rotational period of the neu- 
tron star, or the 35 day intensity variation caused by a 
precessing, inclined, warped, and twisted accretion disk 
(see [Pettcrson 1975| ; |Schandl fc Meyer 1994] ) . 

During the 35 day cycle the X-ray intensity shows 
strong variation, with two maxima in intensity called the 
"main-on" and the "short-on" . This intensity variation is 
caused by periodic covering of the neutron star by the 
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Figure 1. From top to bottom: RXTE PC A, RXTE ASM, 
and CGRO BATSE light curve over the time of the turn- 
on. The X-ray flux is steadily increasing from observa- 
tion 00 until observation 31. During the time of the turn- 
on an eclipse, anomalous and pre-eclipse dip occurred. For 
a further description see text. 



outer (main-on) or the inner edge (short-on) of the accre- 
tion disk. The "turn-on" marks the beginning of the main- 
on and therefore the beginning of the 35 day cycle. During 
the turn-on the outer edge of the accretion disk steadily 
opens the line of sight to the neutron star. While the outer 
accretion disk rim moves out of the line of sight, the ob- 
served X-ray intensity steadily increases, until it reaches 
the maximum after about 2-3 days. Fig. |] shows the X- 
ray light curve of a complete turn-on observed in 1997 
September with RXTE PCA. In this paper we present 
preliminary results of our combined spectral and tempo- 
ral analysis, and resulting simulated pulse variation over 
the turn-on. 



2. Turn-on light curve and spectral analysis 

The complete X-ray light curve observed with RXTE is 
shown in Fig. [j]. During the turn-on the observed flux 
steadily increases until its maximum in orbit 31. Around 
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Figure 2. Evolution of the pulse profile in the energy range 
from 9.0 to 13.0 keV over the time of the turn-on. From 
top to bottom: Pulse profile of orbit 30, 23, 20, 17, and 04- 
The pulse profile of orbit 30 (red) is used as a template for 
the pulse simulations. All profiles are normalized relative 
to the maximum flux in the profile of orbit 30, which is 
set to unity. Pulse phase zero is defined as the maximum 
flux in the pulse profile, thus absolute phase information 
is lost. 

MJD 50704.8 and MJD 50705.3 an anomalous dip, a pre- 
eclipse dip, and at MJD 50705.5 an eclipse occurred. For 
our spectral modeling and pulse analysis we ignored the 
times of the dips (orbit. 10-14) and the time of the eclipse 
(orbit 19-22) because our spectral model is not applicable 
during these times. The numbers in the top panel of Fig. [l] 
identify individual RXTE orb its starting from to 31. 

As presented earlier (e.g., Kuster et al. 1999 ), a spec- 
tral analysis and an analysis of the pulse profile reveals 
that, 

— the spectral evolution during the turn-on can be de- 
scribed by a combination of direct, absorbed, and scat- 
tered radiation, 

— the intrinsic shape of the pulse profile is not changing 
over the time of the turn-on, but the pulse signature 
is steadily washed out towards the beginning of the 
turn-on (see Fig. ||) , and 

— the pulse profile and the spectrum at the end of the 
turn-on (later than orbit 30) shows no significant sig- 
nature of absorbed or scattered radiation. 

These results are in agreement with earlier observations 
of e.g. pavison fc Fabian (1977)|, [Becker ct al. (1977) 



or Parmar et al. (1980), who found the contribution of 
two components to the overall observed X-ray spectrum: 
a "scattered" component and a "absorbed" component. 
The origin for the scattered com ponent is generally in- 
terpreted (e.g. Mihara et al. 1991 ) as the influence of an 
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Figure 3. iV c . a b s and iV ej ; on over the time of the turn- 
on. The astericks mark iV a b s , c determined via the spectral 
analysis using an absorption model. The diamonds mark 
^Ve,ion we get from our pulse analysis. This -/V c ,ion can be 
attributed to a scattering medium which modifies the shape 
of the pulse profile. 



extended, ionized accretion disk corona at a distance of 
r < 5 x 10 11 cm. Recent observations with XMM-Newton's 
RGS show strong indication for the presence of an accre- 
tion disk corona in Her X-l as well (Burwitz et al. 2002). 
The "absorbed" spectral component is attributed to the 
influence of the outer accretion disk rim partially obscur- 
ing the central emission region. 



3. Simulating pulse variation 

Taking the results from the spectral analysis into account, 
we can assume that the variation of the pulse profile at 
the beginning of the turn-on is caused by scattering and 
photoelectric absorption in the medium covering the line 
of sight to the neutron star (the outer accretion disk rim 
and corona). Furthermore we can assume that at the end 
of the turn-on, when the accretion disk rim has moved 
out of the line of sight, the pulse profile represents an 
unaffected emission characteristic. Thus we can use the 
pulse profile observed in, e.g., orbit 30 as a template and 
investigate the effects of a scattering and absorbing corona 
on the pulse shape depending on Nh = N e and the size of 
the scattering region. 

For our simulations we used a Monte Carlo code devel- 
oped at our Institute, assuming simplified slab geometry 
and taking into account Compton scattering, photoelec- 
tric absorption, and fluorescent line emission. As a result 
we obtain the Green's function G(t,to) as a function of 
Nn=N e , the thickness of the slab d, and the energy band, 
and angle-dependent photon spectra of photons leaving 
the slab. We simulated G(t,to) for two cases, a fully ion- 
ized corona and a neutral corona. The total Green's func- 
tion for a partially ionized corona can then be calculated 
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Figure 4- Fit results for orbit 08. The dashed line marks the observed pulse profile and the solid line the simulated pulse 
profile. The simulated pulse profile is a modified version of the pulse profile of orbit 30 (see Fig. ^). 



as Gtot(i) to) — (1~ a ) Gion+a Gabs, where a is the fraction 
of ionization. For G a bs we used A^abs determined by the 
spectral analysis. With a variable A^km for Gi on we can 
then use this simulated Green's function to calculate mod- 
ified pulse profiles an observer located at infinity would see 
according to: 

ri 



I°°(t) 



G{t,t )I(to)dt 



(1) 



where I (to) is the intensity of the template pulse profile 
at time t . Finally we perform a \ 2 minimization fit of 
simulated to observed pulse profiles for each single RXTE 
orbit shown in Fig [l]. As an example, the result of the fit 
for orbit 08 is given in Fig. ^. All four energy bands were 
fitted simultaneously. By combining this procedure with 
a spectral analysis based on a partial covering model (see 
Kuster et al. 1999), we are able to separate the amount 



of scattered and absorbed radiation from purely scattered 
radiation, and determine the thickness d of the corona as 
well as the uncertainties, respectively. 

4. Results 

Preliminary results of our pulse analysis are shown in 
Fig. ||. The crosses in Fig. || show the evolution of neu- 
tral density A^h = A^abs deduced from the spectral anal- 
ysis. As reported earlier by us (Kuster et al. 2001), N, 



abs 



is increasing at the beginning of the turn-on and starts 
to decline orbit 05. This behavior can be understood due 
to an additional flux contribution of a scattering corona 



which comes important at high densities and thus dom- 
inates the spectrum during early phases of the turn-on. 
The diamonds in Fig. [| show the results we get from our 
pulse profile analysis. This A^ion can be attributed to the 
flux contribution of a scattering corona moving slowly out 
of the line of sight to the neutron star. 
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